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Abstract. We have derived the mean broadening profile of the star V 102 in the region of the open cluster IC 4665 from high 
resolution spectroscopy. At a projected equatorial rotation velocity of v sin i = (105 ±12) km s -1 we find strong deviation from 
classical rotation. We discuss several scenarios, the most plausible being strong differential rotation in latitudinal direction. For 
this scenario we find a difference in angular velocity of AC1 = 3.6 ± 0.8radd~' (AQ/fl = 0.42 ± 0.09). From the Ha line we 
derive a spectral type of A9 and support photometric measurements classifying IC 4665 V 102 as a non-member of IC 4665. 
At such early spectral type this is the strongest case of differential rotation observed so far. Together with three similar stars, 
IC 4665 V 102 seems to form a new class of objects that exhibit extreme latitudinal shear in a very shallow convective envelope. 

Key words, stars: rotation - line: profiles - stars: individual(IC4665 V102/P35) - clusters:individual(IC4665) 



1. Introduction 

The substantial difference between photospheres of solar-type 
stars and A-type stars is the existence of a convective enve- 
lope. Late-type stars harbor convective envelopes where tur- 
bulent motions of photospheric plasma can occur. In convec- 
tive envelopes, differential rotation can trigger a magnetic dy- 
namo giving reason for the observed plethora of activity phe- 
nomena. The onset of convection may be traced by the onset 
of co r onal X-ray emission a round spectral type A7 (Schmitt, 
ll997l) . lG7a v&Nagel ( 1989) searched for the onset of convec- 
tion analyzing line bisectors of slowly rotating stars. In their 
targets a bisector reversal was found around spectral type F0. 
Stronger asymmetries were found in the stars at the hot side of 
the boundary indicating higher photospheric velocities. 

While stellar rotation velocities have been measured for al- 
most a century, differe ntial rotation is more difficult to detect 
fe.g. lGravul977l[l998l and references therein). Measurements 
of differential rotation have been reported in stars rotating 
faster than v sin/ > lOkms 1 a nd covering all spectral types 
harboring convective en velopes jReiners & Schmitl l2003albl: 
iReiners & Rovetll2004 . 



* Based on observations carried out at the European Southern 
Observatory, Paranal, 7 1.D-0 127(A) 



2. Data 

Our data of V 102 have been taken on April 18th 2003 with the 
FLAMES multi-object facility at ESO/VLT. Using the UVES 
spectrograph at a resolution of 38 000 we gained a signal-to- 
noise ratio of 150 in a one hour exposure. The spectra were re- 
duced using the UVES context embedded in the ESO-MIDAS 
package. The observed spectral range covers the region from 
4800 A to 6800 A with a gap at 5800 A. In a second observa- 
tion during the same night we collected a spectrum with the 
FLAMES facility using the GIRAFFE spectrograph which cov- 
ers the Can H+K resonance lines. A small region of the UVES 
spectrum is plotted in Fig.^(black spectrum). 



3. V102 and its membership to IC 4665 

IC46 65 is a young (50-70 Mvr: IProsseAll993llGiampapa et all 
1998) open cluster at a distance of roughly 350 pc. The 
first photometric meas urements of V 102 were done by 
iMcCarthv & 0'Sullivanllll969l) . Although V 102 fitted onto the 
main sequence in their color-magnitude diagram, they argued 
that V 102 is a "less probable" member of IC 4665 due to their 
proper motion measurements. More than twenty years later, 
Prosser (1993) carried out a photometric study of IC 4665 and 
derived more accurate proper motions. From that he derived 



2 



A. Reiners et al.: Strong latitudinal shear in IC4665 V102 



1.00 



X 

E 0.95 

XI 
01 
N 

13 

a 

O 

c 

0.90 



0.85 




5550 



5650 



5600 
Wavelength [A] 

Fig. 1. Data (black) and fit (grey) in a small part of the used wavelength region. In the inset the two broadening functions derived 
from the spectral regions 4950 A - 5550 A (grey) and 5440 A - 5745 A (black) are shown. Both broadening functions perfectly 
match one another. 



an 86% membership probability for V 102. However, he also 
carried out CCD-photometry and found V 102 not to fit onto 
the cluster main sequence according to its brightness in V. The 
results from photometric measuremen ts in V 102 are given in 
Table d The V -magnitudes de rived bvlMcCarth v & O'Sullivanl 
( 1969) and by Prosser ( 1993) are significantly different while 
B-magnitudes are iden tical and are al so reproduced by the 
value measured bv lKislvuketaD (ll999). Since magnitudes in 
the B-filter are identical in all three references we suspect that 
one of the magnitudes in V is incorrect, giving rise to the dif- 
ferent values of B - V. 

The spectral types according to the B - V measurements are 
either A9 or F5. In this spectral range the Ha-line is very sen- 
sitive to temperature and we can deduce the spectral type from 
our data. The Har-line of V 102 is plotted in Fig.|2]together with 
spectra of HD 115 810 1 (B - V - 0.27, v sin/ = 100 kins' 1 , 
light grey) and HD 120987 2 (B-V = 0.44,v sin/ = 8.5 kms 1 , 
dark grey) spun up to v sin / m lOOkms -1 . From the excellent 
fit to the Ha-line of HD 1 15 810, we conclude that the color of 

V 102 must be similar to HD 1 15 810. W e thus adopt the value 
of B - V = 0.26 given in Prosser ( 1993), i.e., from photometry 

V 102 is likely t o be a background star and not a member of 
IC 4665 dProsseilfl993l) . 



Table 1. Photometric data of IC 4665 V 102. 



V 


B 


B-V 


V-I 


B-R 


designation 


11.65 


12.09 


0.44 






V 102 a 


11.82 


12.08 1 


0.26 


0.41 




P35 b 




12.06 






0.69 


354 699° 



Calculated fro m V and B-V, n ot given in the catalog 
McCarthy & Q'SullivanHl969l) 



ProsseJll993l 
iKislvuk et alj I 



1 taken with the FOCES spectrograph, DSAZ, Calar Alto 

2 taken with the FEROS spectrograph, ESO, La Silla 



<1999li 



4. Broadening Function 

From our observations we derived the me an line br oadening 
profile following the strategy outlined in Rei ners & Sch mitt 
(l2002al) . Normalization of the spectrum was done by eye, even 
at rapid rotation the continuum is still visible. Using absorp- 
tion line wav elengths from the Vienna Atomic Line Database 
(Kun kaet all fl 999) as a reference, we iteratively determine 
the shape of the broadening function and equivalent widths of 
all incorporated lines. To check against systematic errors we 
performed the fit in two overlapping wavelength regions. The 
wavelength region at 4950 A - 5550 A incorporates 185 lines, 
the second region at 5440 A - 5745 A 137 lines. 
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Table 2. Results derived from the broadening function 
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Fig. 2. Ha-line of V 102 (black). Overplotted are spectra of 
HD 115 810 (B-V = 0.27; light grey) and HD 120 987 (B-V = 
0.44 spun up to v sin i * 100 km s -1 ; dark grey). The spectrum 
of HD 115 810 resembles the one of V 102 much better than 
the spectrum of HD 120987 does. 



We calculate the variances of the broadening function by 
taking into account the uncertainties of the two nearest neigh- 
bors of each pixel. We overplot in Fig.fl]our final fit to the data 
in grey. The line shapes of all lines are well resembled within 
data quality. In the inset of Fig.^we compare the two broaden- 
ing functions independently derived from the two overlapping 
wavelength regions mentioned above. Both broadening func- 
tions resemble each other indicating indepence on the wave- 
length region; the broadening function with its errors is plotted 
in black in the left panel of Fig. [3] Although every pixel was a 
free parameter during the fit, the broadening profile is symmet- 
ric within the uncertainties. 

Quantitative conclusions can be drawn from the Fourier 
transform of the broadening profile which is shown with its er- 
rors in the right panel of Fig. [3] In Fig. [5] we also plot the broad- 
ening profile of HD 46273, a star of similar spectral type and 
v sin / (F2V, v sin/ = 107 km s -1 ). This profile was achieved 
applying the same procedure to data obtained with the FEROS 
spectrograph at ESO/La Silla dReiners & Schmittl l2003bl) . it 
has a shape as expected from classical solid body rotation. 



?1 



a = AQ/Q. 



An 



(105 ± 12) km s" 1 1.39 + 0.01 0.42 ± 0.09 3.6 ± 0.8radcT 



5.1. Differential rotation 

Convenient quantities characterizing stellar rotational broaden- 
ing are the first two z eros q\ and q? of the broadening profile's 
Fourier transform (cf. lReiners & S chmitt. 2002a). The value of 
v sin / can be derived from q\. Furthermore, in a rigid rotator, 
the ratio q2lq\ always is larger than 1.72. For V 102 we derive a 
ratio qilqi — (1.39 ±0.01) < 1.72. This small ratio is visible in 
the narrow first side-lobe of the Fourier transform compared to 
the profile of HD 46 273 in the right panel of Fig. [3] To quantify 
differential rotation we assume a rotation law parameterized in 
analogy to the solar case: 



Q(/) = Q E quator(l ~ Of SHI /), 



(1) 



with / heliographic latitude and Qriquator angular rotation speed 
at the equator. Differential rotation can be described by the pa- 
rameter a, which is the difference of the equatorial and polar 
rotation velociti es relative to the equatorial v elocity. 

As shown in Rein ers & Schmittl J2003al) . this parameter a 
can be calculated from the observed ratio qijq\, specifically 
qilqx K oil Vsin /, with i the inclination angle. For the case 
of V 102 we derive a = 0.42 ± 0.09, and with a radius of 
R - 1.5 Rq the difference between polar and equatorial angu- 
lar velocities is A = 3.6 ± 0.8radd _1 . 

Such strong a shear in a fast rotator is expected to maintain 
a strong magnetic dy namo and we would ex pect significant X- 
ray flux from V 102. Giampapa et al. ( 1998) carried out X-ray 
observations of IC 4665 with no detection of X-ray flux from 
V 102. However, for a distance of d — 350pc, we estimate 
an upper limit of X-ray luminosity Lx S= 3.4 10 29 erg s _1 . This 
is abo ut the typical X-ray lu minosity detected in that spectral 
range (Pallavi cini et all Il98ll) . Since V 102 is likely to be a 
background object, the X-ray observations do not contradict 
the picture of a strong magnetic dynamo being driven by rapid 
differential rotation. 



5. Interpretation 

The broadening function of V 102 shows significant differ- 
ences to classical rotation broadening. Within the error bars 
a difference is clearly visible between the profiles of V 102 
and HD 46273. At a projected equatorial rotation velocity of 
v sin/ = 105 km s -1 , the broadening functions are clearly ro- 
tation dominated. Since the derived broadening functions are 
symmetric, it is unlikely that they are dominated by starspots 
and we argue that the profiles are due to the rotation law of the 
stellar surface. In the following, we will discuss the amount of 
differential rotation we derive for V 102 and the consequences 
for the surface shear. However, we will also discuss several 
other scenarios that could be reason for such a profile. 



5.2. Alternative scenarios 

Although we believe that differential rotation is the most prob- 
able interpretation of the derived broadening profile, other 
mechanisms can be potential reasons for the deviations from 
classical rotational broadening as well. We discuss three alter- 
native scenarios. 

5.2.1 . Low inclination 

Very rapidly rotating stars become gravity darkened, i.e., they 
become cooler at the equator. At (unprojected) rotational ve- 
locities larger than 150km s _1 , the shape of the broadening 
profile, and espec ially the ratio qi/qi, is affected by grav- 
ity darkening (cf. iReinersl 120031) . From the models given in 
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Fig. 3. Broadening functions of IC 4665 V 102 (black) and HD 46 273 (grey) in velocity (left panel) and Fourier space (right 
panel). The profile of HD 46273 is resembles classical rotation broadening. Significant differences are visible to the profile of 
IC 4665 V 102. Note especially the narrowed first sidelobe in the Fourier transform - a typical signature for differential rotation 
with the equator rotating faster then the pole. 



lReinersN2003l) we can calculate the equatorial velocity needed 
for our observed ratio of qilqi = 1.39; we get a value of 
v ro t = 800kms _I , which is far beyond breakup velocity. Thus 
we can rule out gravity darkening as the main reason for the 
observed profile. 

5.2.2. Bright Polar Spot 

The surface of V 102 may be covered with spots. If the star 
were magnetically active, one would expect cool spots cover- 
ing the star. Cool spots generate characteristic signatures in ab- 
sorption profiles. Since the derived broadening profile is sym- 
metric, the surface coverage of spots must also be symmetric 
to a high degree. One configuration that would satisfy these re- 
quirements is a polar spot. We examined the influe nce of spots 
on the broadening profile (cf. Reiners & Schmitt, 2002b) and 
found that cool polar spots do change the profile shape and the 
ratio qilq\ but in the opposite direction; qilqx becomes larger. 
A profile shape as the one observed could only be due to a 
hot polar spot. From our models we estimate that a polar cap 
with a radius larger than 30° at a temperature about 1000K 
hotter than the rest of the star is needed to produce such a pro- 
file shape. This could be similar to what is observed in T Tauri 
stars accreting material at the poles, and one obvious sign of 
such accretion would be Ha emission. In V 102, no sign of Ha 
emission can be found (cp. Fig.[5J- 

We emphasize that this may nevertheless be a scenario al- 
ternative to differential rotation, but we see no physical connec- 
tion to any phenomenon observed on other stars. 

5.2.3. Composite spectrum 

The derived broadening profile could tentatively also be a com- 
position of two spectra from rigidly rotating stars. Especially 
the curvature at around ±70kms _1 could be interpreted as an 
indication for that scenario (although it can also be explained 
by differential rotation). The profile could then be interpreted 
by a spectrum of the first star rotating at v sin / « llOkms -1 , 
superimposed by a spectrum of a second object rotating at 
v sin i * 70kms _1 . Thus both stars would have either sig- 



nificantly different rotation velocities or were seen under dif- 
ferent angles of inclination, i.e., a hypothetical double sys- 
tem could not be corotating. Since we derived identical pro- 
files from two different wavelength regions covering several 
hundred Angstroms, both stars had to be of similar spectral 
type, otherwise the line profile would change at different wave- 
lengths. In order to resemble the derived profile, both stars 
must also contribute approximately the same amount of flux. 
Given similar spectral types of the components, that implies 
that both stars are at the same distance and are probably physi- 
cally bound. The double system would then emit roughly twice 
the flux a single star would, i.e., the real magnitude of each 
member would be 0.75 mag higher than the measured ones 
putting V 102 even further away from the main sequence in 
the color-magnitude diagram. 

In this scenario, V 102 would be a physically bound system 
consisting of two stars of similar spectral type and different ro- 
tation velocities located far outside the cluster but comoving 
with the cluster. We would expect such a spectrum to be vari- 
able. Although we consider this scenario very unlikely, we aim 
to reobserve the star to check the profile for variability. 

6. Conclusions 

We derived the mean line broadening profile of V 102, a back- 
ground star of spectral type A9 in the region of the young open 
cluster IC 4665. From the broadening profile we determined a 
projected rotation velocity of v sin; = 105 ± 12kms _I . The 
broadening profile significantly differs from classical rotation 
broadening and we discuss several scenarios for such a pro- 
file. The most plausible scenario is that V 102 does not rotate 
as a solid body but that its equator rotates at a higher velocity 
than the pole. For that case we derive a difference in polar and 
equatorial angular velocities of AO = 3.6 + 0.8radd~', or 
AQ/Q. = 0.42 ± 0.09. 

This is so far the strongest observed case of differential 
rotation. The equator of V 102 laps the polar region once 
every 40 hours, compared to roughly 120 days on the Sun. 
Three other stars with strong differential rotati on of the orde r 
of AQ/D. = 0.3 have been reported in Reiners & Rover 
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they have colors of B - V — 0.29 or 0.30. Our result supports 
evidence for differential rotation being strongest in a class of 
rapidly rotating objects with very shallow convection zones. 
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